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Abstract 

Assume that f(s) = F'(s) where F is a double- well potential. Un- 
der certain conditions on the Lipschitz constant of / on [— 1 , 1] , we 
prove that arbitrary bounded global solutions of the semilinear equa- 
tion Au = /(it) on hyperbolic space H™ must reduce to functions 
of one variable provided they admit asymptotic boundary values on 
gn-i _ ^jjn which are invariant under a cohomogeneity one sub- 
group of the group of isometries of EP. We also prove existence of 
these one-dimensional solutions. 



1 Introduction 

The equation —Au + u(u 2 — 1) = is the Euler-Lagrange equation for the 
functional f ^\Vu\ 2 +^(u 2 — l) 2 ; it has been intensively studied for many 
years, and has many applications in geometry, physics and elsewhere. The 
important feature of this energy functional for many of these investigations 
is the 'double- well' behaviour of the nonlinear term F(u) = (1/4) (u 2 — l) 2 , 
and indeed, many results generalize directly to extremals of other energy 
functionals where F is replaced by a more general function with qualitatively 
similar behaviour. 

One important theme in this research is to study the symmetry prop- 
erties of solutions, in particular those which are defined globally on R n . 
Most famous here is the De Giorgi conjecture, which asserts that, so long as 
n < 8, any solution u(x) of this PDE which takes values in [—1,1], and sat- 
isfies certain monotonicity properties to be made explicit below, necessarily 
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has the form U(£ ■ z), where i £ M. n is constant, z £ M. n and U(s) satisfies 
U" = F'(U); in other words, u reduces to a function of one variable. This 
has been proved for n < 3 [U |2J. The cases 4 < n < 8 have been proved 
by Savin [12], but with an additional assumption beyond monotonicity that 
the solution converges to ±1 at infinity. For related results see [HElSj, and 
[7] contains a survey. 

Our goal here is to investigate these symmetry results from a broader 
perspective. Let (X,g) be a complete homogeneous Riemannian manifold, 
and G = Isom (X) its group of isometries (which thus acts transitively) . 
Consider bounded solutions u of the equation 

A g u - F'(u) = 0, (1.1) 

where F is a double-well potential function. More precisely, 

F G C 2 (R), F > 0, F(s) -> oo when \s\ -► oo, , . 

F- x (0) = {-1, +1} and F"(±l) > 0. 1 j 

We shall also always assume that sF'(s) > when |s| > 1, which is a 
technical convenience which (by the maximum principle) allows us to assume 
only that u is bounded rather than that u(z) £ (—1, 1) for all z; furthermore, 
we also assume that the only other critical point of F is a local maximum 
at s = 0, though this too can be weakened. The main question is as follows: 

If u is a bounded solution of (jl.ip defined on all of X, then under 
what conditions is u necessarily invariant under some subgroup 
of positive dimension of G? In particular, under what conditions 
is u invariant under a subgroup H C G which has orbits of 
codimension one in X? 

The constant solutions u = ±1 are invariant under G itself, but we are more 
interested in solutions which 'jump states', i.e. such that u — ► —1 in certain 
directions at infinity and u — ► +1 in others. 

The types of hypotheses we shall place on u concern its 'asymptotic 
boundary values', i.e. its limits in different directions at infinity. This differs 
from the hypotheses in the De Giorgi conjecture where one assumes only 
that u is monotone along some fixed family of parallel lines; our hypothesis 
is analogous to the so called Gibbons condition [3]. Part of the work below 
involves showing that once one has imposed certain asymptotic boundary 
values, the solution necessarily has a monotonicity property; however, as we 
indicate below, for certain choices of boundary values, if F"(0) is sufficiently 
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large, there are ODE solutions of (jl.ip which are not monotone but do have 
these correct asymptotic boundary values. 

We shall not formulate the precise problems we study in this generality, 
but consider only the particular case that X is the n-dimensional hyperbolic 
space HP, and hence G = Isom (HP) = 0(n, 1). This is of interest since the 
isometry group G is substantially more complicated than the group E(n) of 
Euclidean motions. Many of the results here can be proved in much greater 
generality; for example, Theorem 2 below remains valid for any rank one 
symmetric space X (and could even be formulated suitably for higher rank 
symmetric spaces of noncompact type). The key property needed is that G 
contains a sufficiently large collection of reflections. However, for simplicity 
we focus on this problem on hyperbolic space. 

The geodesic compactification HP = HP U 5 n_1 is obtained by adding a 
point at the end of each geodesic ray emanating from some fixed point. This 
asymptotic boundary has a natural topology and conformal structure, and 
isometries of HP induce conformal transformations of the sphere at infinity 
and conversely. A function u on HP has asymptotic boundary value eft on 
S"™ -1 if u is continuous (a.e.) on HP and its restriction to the boundary is 
equal to (p. 

Briefly, we consider solutions u of (11. ip on HI™ with asymptotic boundary 
values invariant with respect to a cohomogeneity one subgroup H C G, and 
the main question we ask is whether the solution u itself is also invariant 
with respect to H. Since HP/if is one-dimensional, this means that u is 
determined by a function of one variable (which hence satisfies an ODE). 

Before stating our results precisely, we recall the three main models for 
HP: the upper-half-space, 

Rl = { Z = ( x ,y) G R n : x > 0, y G BP -1 }, g = ^ + ^ 
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the Poincare ball, 

B n = {zeM. n :\z\<l}, g= (1 4 ^ 2 2)2 , 

and via geodesic normal coordinates, 

1" = {z = rto : r > 0, lo G S 11 ' 1 }, g = dr 2 + sinh 2 r duJ 2 . 
There are three subgroups of G which have orbits in HP of codimension one: 
fle = 0(n), H p = E(n), H h = 0(n-l,l). 
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The subscripts 'e', 'p' and 'h' stand for elliptic, parabolic and hyperbolic, 
respectively. The action of H e is by rotations around a fixed point p € HP 
or by reflections across a totally geodesic hyperplane containing that point; 
in the ball model, with p at the origin, or in geodesic polar coordinates, this 
corresponds to the standard rotation group acting on B n or HP, respectively. 
Next, Hp acts by Euclidean motions on each horosphere centered at a point 
q £ S 11 ^ 1 = <9HP; in the upper-half-space model, with q the point at infinity, 
these horospheres are the hyperplanes x = constant and the action appears 
as the standard linear one in each HP -1 . Finally, the inclusion of in 
G depends on the choice of a totally geodesic hyperplane P = HP -1 C 
HP. Any isometry h of P extends (uniquely up to a reflection across P) 
to an isometry of the entire space HP, and the subgroup H^ consists of all 
isometries obtained in this way. In the ball model, we could take P = {z n = 
0}, while in the upper half-space model a convenient choice is P = {\(x, y)\ = 
1}. The hypersphere P n <9HP = S n ~ 2 C S n ~ 1 divides the boundary sphere 
into two components, which we denote by S- and S+, respectively. Finally, 
we let t be the function on HP which gives the signed distance from P, so 
t > on one component of HP \ P and t < on the other. 

Any function on S 71-1 = <9HP invariant under the induced action of H e 
must be constant; a function invariant with respect to the action of H p 
can assume only two values, one at the parabolic fixed point q and the 
other on the other orbit, S 11 ^ 1 \ {<?}; finally, any function on S 11 ^ 1 which 
is invariant under the subgroup of H^ which does not interchange the two 
components of HP \ P takes on one constant value on all of S- and another 
on all of S+. We present two main theorems. The first guarantees the 
existence of Hi invariant solutions for each of these subgroups, i = e,p,h, 
and the second gives conditions under which all solutions (with i?j-invariant 
boundary values) must themselves be Hi invariant. This second theorem 
includes some decay hypotheses on the solution, and to show that these are 
necessary we also prove the existence of non-symmetric solutions which have 
just slightly weaker decay. 

Before stating these, we introduce a few important constants associated 
to this problem. Let L denote the Lipschitz constant of / on the interval 
[—1, 1], and A = — /'(0) > 0. Consider the two equations: 

a 2 - (n- l)a + A = 0, a 2 - (n - l)a + L = 0, (1.3) 

and denote the roots of these by a± and f3±, respectively. Note that it is 
always true that A < L. We shall need to assume at various points below 
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that L<(n- l) 2 



/4; this ensures that all of these roots are real, and satisfy 



< a_ < (3_ < 



n — 1 



< f3 + < a + < n - 1. 



(1.4) 
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Theorem 1.1 Consider solutions of which are invariant under H e , 

Hp or Hh, where f = F' and F satisfies the hypotheses in Itl. 

i) When i = e, the only ODE (i.e. radial) solution which has constant 
asymptotic boundary value ±1 is the constant solution u = ±1. There 
is a one-parameter family of radial solutions, all of which decay like 
(1 — as \z\ — * 1 in the ball model; however, if \u\ < C(l — \z\) a 

for some a > a_, then u = 0. 

ii) When i = p, there exists an Hp-invariant solution u p which is unique 
up to translation, and which tends to on \ {q} and to +1 at 

{q}, where q is the parabolic fixed point of H p . For generic functions 
f, this solution decays like x a ~ as x \ 0. On the other hand, there 
exists no Hp-invariant solution which tends to —1 along M n_1 and +1 
at infinity, or vice versa. 

Hi) When i = h, there exists an H^-invariant solution which is mono- 
tone as a function of the signed distance t, and which tends to —1 in 
S- and +1 in S + . 

Note that we are not asserting the uniqueness of the ODE solution in case 
iii), and indeed it is likely that there is a sequence of 'multibump' solutions 
of this problem which have an arbitrary finite number of oscillations. 
The converse statement is provided by the 

Theorem 1.2 Let u be an arbitrary bounded solution of where f = F' 

and F satisfies the hypotheses in 

i) Let rj = — 1 or +1. // \u(z) — rj\ — > as z — > oo in W 1 , then u = 77. 
On the other hand, suppose that f is Lipschitz with Lipschitz constant 
L < (n — l) 2 /4 and f(s) G C l,s in some neighbourhood of s = for 
some 5 > 0. // |-u(2:)|/(l — \z\) a < C for some a > a_, where z is the 
coordinate in the ball model), then u = 0. It suffices to require only 
that f is Lipschitz if this exponent of decay satisfies a > (3-. 

ii) Suppose that, as before, f{t) is Lipschitz with Lipschitz constant L < 
(n — l) 2 /4. // |u| < Cx a on all of W 1 for some a > a_ (where we 
use the upper half-space variable x), and u — ► +1 in any conic region 
x > C\\y — ?/o I + C2 as (%, y) — > 00, then u is invariant under H p , and 
hence is a function of x alone. 
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Hi) Ifu — > — 1 uniformly on compact subsets of S- andu — > +1 on compact 
subsets of S+, then u depends only on the signed distance t from the 
totally geodesic HP 1-1 . 

If u is an arbitrary solution which is bounded by x a for a > a_, then 
by part ii), u is invariant under H p ; however, by part ii) of the preceding 
theorem, this ODE solution does not decay so quickly as x \ 0, at least 
for most choices of function /. The conclusion is simply that there are no 
solutions with this rate of decay. 

As was noted earlier, unlike for the De Giorgi conjecture in W 1 , we 
assume conditions about the asymptotic boundary values but do not impose 
monotonicity conditions on the solution. 

Our final result is to show that the decay hypotheses in parts i) and ii) 
of Theorem 1.2 are necessary. 

Theorem 1.3 Let f be as above. 

i) There is an infinite dimensional family of solutions of §1.1}) . each el- 
ement of which decays like (1 — \z\) a ~ as \z\ — > 1 in the ball model of 
HP. 

ii) There is an infinite dimensional family of solutions of hl.l\) . each el- 
ement of which is bounded by x a ~ in the upper half-space model of 

The equation (jl.ip has many other interesting global bounded solutions 
on EP. It would be interesting to construct bounded solutions which are 
piecewise constant on the sphere infinity, taking only the values 0, ±1 there. 
The results here address only the very simplest situations of this sort. It 
would also be interesting to understand if monotonicity but no assumptions 
on the asymptotic boundary conditions are enough to imply that the solution 
is one-dimensional. We shall come back to these issues elsewhere. 

2 ODE solutions 

We now prove the existence of one-dimensional solutions to (II. ip which are 
invariant under either H p or H^, and discuss the uniqueness of solutions in 
the H e case. 

2.1 The elliptic case 

The existence of radial H e invariant solutions is obvious, since we can simply 
let u = or ±1. We now consider the uniqueness of these ODE solutions. 
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Proposition 2.1 The only H e solution of fZHP which tends to ±1 at in- 
finity is the constant solution. However, there is a one-parameter family of 
radial solutions to this equation which tend to at infinity. 

Proof: For radial functions and in geodesic normal coordinates, (jl.lj) takes 
the form 

u + (n — 1) coth r u = f(u). (2-5) 

Suppose that u(r) — > ±1 as r — > oo. Multiplying by v! and integrating from 
to A yields 

1 f A 

-u'(A) 2 + (n - 1) / coth ru'(r) 2 dr 

2 Jo 

A 

f(u(r))u'(r) dr = F(u(A)) - F(u(0)), 

where we used that u'(0) =0. As A -> oo, u'(A) -> and -> 0, 

whereas F(u(0)) > 0, so taking the limit gives 

POO 

(n-1)/ coth r u'(rf dr = -F(u(0)) < 0; 

this implies that u' = 0, hence u = ±1. 

This calculation does not provide any information when u — » at infinity 
since 1^^(11(74)) > 0. In fact, there is a one-parameter family of noncon- 
stant radial solutions tending to as r — > 00. To see this, consider the solu- 
tion u a (r) of (|2.5h which satisfies the initial conditions u o (0) = 06 (-1, 1), 
u' a (0) = 0. By the maximum principle, |it a ( r )| < 1 for all values of r for 
which the solution exists, and because of this boundedness it is easy to see 
that u' a and u" a also remain bounded, so that u a (r) exists for all r G M + . 

We next show that u a has a limit as r — > 00. If it were monotone 
there would be nothing to prove otherwise there would exist a sequence 
{rj} tending to infinity where u' a vanishes. There would in fact exist some 
subsequence r'- such that u a (r'j) tends to a limit (and u' a (r'-) = for all j). 
Multiply (|2.5p by u' a and perform exactly the same integration by parts as 
above on the interval (r^-,r^ +1 ) to get 

(n - 1) / 1+1 coth ru' a {rf dr = F{u a {r' j+1 )) - F(u (r^)) -» 0. 

Since the left hand side is nonnegative, the sequence {F{u a {r'-))} is mono- 
tone increasing. If the lengths of all the intervals — r'-\ — > 0, then 
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since u' a is bounded, we conclude that u a itself has a limit. Otherwise, we 
conclude at least that the L 2 norm of u' a is uniformly small on each such 
interval, and hence it converges uniformly to too, here too we conclude 
that u a converges. 

Finally to prove that u a converges to zero, writing (|2.5p 

[g(r)u' a ]' = g(r)f(u a ) 

with g(r) = (sinhr) n_1 and integrating it on each (r^-,r^ +1 ) gives 

0= I 3+1 g(r)f{u a {r))dr. 



This implies that for each j there exists r OJ - £ (r'p r'j +1 ) such that f(u a (r j)) = 
0. This gives that u a converges to either or ±1, and since we already 
proved that it cannot converge to ±1 without being identically equal to ±1, 
we conclude that u a (r) — > as r — > oo, as claimed. □ 

2.2 The parabolic case 

Consider H n in the upper half-space model, placing the parabolic fixed point 
at infinity. An fZp-invariant solution satisfies the equation 

x 2 u"(x) + (2 - n)xu'{x) - f(u) = 0, (2.6) 



where x £ IR . We assume that / satisfies the conditions in (II. 2p . and 
moreover that /'(0) < 0. 

Proposition 2.2 There are no bounded solutions of \2. 6\) which satisfy 
lim u(x) = —1, lim u(x) = 1. 

x— >0 a^+oo 

On the other hand there exists a solution, unique up to the rescaling x i— > cx, 
c > 0, such that 

lim u(x) = 0, lim u(x) = 1. 

x^0 x^+oo 



Proof: Setting £ = logx, and letting primes now refer to derivatives with 
respect to £, then (12. 6p becomes 

u"(0-(n-lV(0-/(«)=0. (2.7) 
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Setting u' = v yields the equivalent autonomous system 



u' = v 



(2.8) 



v' = 



(n -!)« + /(«) 



which we shall analyze by phase plane methods. There are three critical 
points: (0, 0), which is either an unstable node or an unstable spiral depend- 
ing on whether — /'(0) is less than or equal to, or greater than, (n — l) 2 /4, 
and (±1,0), which are both saddle points with unstable eigenvectors of the 
form (1, a±), where both a± > 0. We are thus looking for heteroclinic orbits 
connecting any two of these critical points. 

Let us first examine the integral curve which emanates from (—1,0) 
as £ — > — oo and moves into the upper half-plane v > 0. There are only 
three possible scenarios: either this orbit crosses the line v = in the open 
interval where < u < 1, or it tends to (1,0) as £ — > +oo, or finally that 
it passes above this point. Note that the orbit cannot cross in the interval 

— 1 < u < since the vector field points upward there. We shall show that 
the first two possibilities cannot occur, which then proves the first assertion 
of the theorem, as well as the second since it then follows that precisely one 
of the integral curves emanating from (0,0) as £ — * — oo converges along 
the stable direction to (1,0), and this is the solution we want. Note that as 

— f'(0) increases, this solution may spiral around the origin more and more 
times; indeed, once — /'(0) > (n — l) 2 /4, it spirals around infinitely often. 
There is a critical value 7 so that if — /'(0) < 7, then the corresponding 
solution stays strictly positive, but if — /'(0) > 7, then the solution becomes 
negative; hence one should not expect solutions to be monotone unless |/'(0)| 
is sufficiently small. 

Thus suppose that (u(£), ?;(£)) is a solution curve connecting (—1,0) to 
(+1,0). Multiply (j2T!) by v! and integrate from -A to +A to get 



since F(—l) = F(+l). This proves that such a solution cannot exist. If 
the integral curve emanating from (—1,0) reaches the n-axis in the range 





F(u(-A)). 



Taking the limit as A — ► 00 yields that 
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< u < 1 at some finite time T, then since u'(T) = v(T) = 0, the same 
calculation yields that 

-(n - 1) f {u'f di = F{u{T)) - F(-l) > 0, 

J — oo 

which is again impossible. The same reasoning shows that there cannot 
be a homoclinic orbit connecting (1,0) to itself. All of this establishes the 
existence of a unique integral curve connecting (0,0) to (1,0). 

Unless / is very special, this solution will not approach (0, 0) along the 
eigenvector with eigenvalue a+, hence for generic /, u(x) ~ cx a ~ as x \ 0. 

The case of a solution going from (0,0) to (—1,0) is handled similarly. 

□ 

There is one specific case where this solution can be written explicitly. 
Let /(it) = ku(l — u 2 ), where k = 2<yn ^ , and set a = ^ip; then 



is an explicit solution of (|2.6jl satisfying 

lim u = 0, lim u = 1. 

x^0 z-+oo 

(Starting the 'time' variable £ at different points along this solution yields 
the family of solutions x a /(b + x a ) for any b > 0.) 

Corollary 2.3 If u is any bounded solution of \l. 1\) such that, in the upper 
half-space model, 

lim u = —1, lim u = 1, 

x^0 x— >oo 

i/ien it does no£ reduce to a function of only one variable. 
2.3 The hyperbolic case 

The solutions of (II. lj) which are invariant under the subgroup -ff^ satisfy the 
ODE 

f7"(t) + (n-l)tanht£/' = /(£/), lim *7(i) = ±1. (2.9) 

t^±oo 

We now prove the existence of such a solution. 

Proposition 2.4 Suppose that f(s) = F'(s), where F satisfies the same 
assumptions as before. Then there exists a solution U of $2.9\) . 
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Proof. We first find a function Ut{P) denned on the interval [— T, T] for 
any T > which minimizes the functional 



over the class of functions u £ H l {[—T, T]; (cosht) n 1 dt) satisfying u(±T) = 
dbl. Since F is nonnegative, it is standard that such a function exists. 
Note that U T (t) satisfies the ODE in fl2J)D and U T {t) £ (-1, 1) for \t\ < 
T. It is also not hard to see that Ut is monotone increasing. Indeed, by 
the maximum principle, no solution of this ODE can have a positive local 
minimum or negative local maximum. Therefore, if Ut is not monotone, 
then at least it must be monotone increasing on two intervals [— T, A] and 
[B, T], where U T < on [-T, A], U T > on [B, T] and U T (A) = U T (B) = 0. 
If monotonicity were to fail, it would be because Ut oscillates on [A, B] . 
However, this cannot be the case for a minimizer, since the function which 
agrees with Ut on [-T, A] U [B, T] and is identically on [A, B] has less 
energy. (There are quite likely nonminimal critical points of this functional 
which do oscillate.) 

In order to prove the existence of hin^^oo Ut, it suffices to show that the 
value Ut(0) remains bounded away from —1 and +1; alternately, we must 
show that the value a = (it where Ut(cit) = remains bounded as T —* oo. 
If this were not the case, consider the sequence of functions Ut(- — clt)', this 
has a limit as T — > oo which is defined either on all of M or else on some 
half-line [—A, oo) or (— oo, B], and has the boundary values or limits —1 and 
+1 at the left and right ends in each case. To be concrete, suppose that 
a,T — * — oo, so that the limit function v(t) satisfies v" — (n — l)v' — f(v) = 
on either 1R or [—A, oo). Multiply this equation by v' and integrate over the 
domain of definition to get 



(the first term on the left is absent if A = oo), which is a contradiction. 
This proves that the limit U of Ut exists, satisfies the ODE (|2.9j) and lies 
in (-1,1) for all t. 

It remains to prove that it has the correct asymptotic limits. Since 
U{t) is monotone increasing, there is some limiting value u = lim^oo U(t). 
However, the eventuality lo < 1 is ruled out by the ODE again since the 
limit of the left side would vanish while lim^oo f(U(t)) = f(u) ^ 0. The 
corresponding statement is true for the limit of U as t — > — oo. □ 



E{u) = [ (hu'(t)\ 2 + F(u(t)) ) (cosh ty^dt 



v'{-A)\ 2 - (n - 1) / (v' f dt = F{1) - F(-l) = 
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We remark that the key property of / required in this proof is that 
f^f(s)ds = 0. 

3 One-dimensional symmetry of solutions 

We now turn to the problem of determining when global bounded solutions 
of (jl.ip are invariant under one of the subgroups H e , H p or H^. The first two 
cases will require an additional decay hypothesis which we shall demonstrate 
later is necessary. This can be stated in several equivalent ways. The most 
convenient is in terms of two different families of generalized eigenfunctions 
of the Laplacian: the plane wave solutions o~ a and the spherical functions 
<p a , which we now define. In upper half-space coordinates, the plane wave 
solution with pole at infinity and exponent a G C is defined by 

<r a (x,y) = x a ; (3.10) 

it satisfies Ao~ a = a(a — (n — l))a a . Thus, when Re a > 0, a a — > when 
x \ and a — > oo when x /* oo. Similarly, for any a £ C there exists 
a unique radial (ff e -invariant) function <j) a on H n which satisfies this same 
eigenfunction equation; when a < ^(n — 1), 4> a ~ (1 — \ z \) a as \ z \ ~^ 1- 

In the ball model, the distance to the boundary sphere in the Euclidean 
metric is comparable to the exponential of the negative of the distance func- 
tion for the hyperbolic metric: 

Ci(i - |z|) < e - rf (^°) < C 2 (l - |^|), 

so it is not hard to express the decay of a a and <j) a in terms of this geodesic 
distance. 

We first present two simple lemmas which will be used in several places 
below. 

Lemma 3.1 Let u be a bounded solution of ( EUP on W 1 . Suppose that Tj 
is any sequence of Mobius transformations. Then the functions u o Tj are 
also solutions of U.l\) and some subsequence of these converges ( at least in 
C 2 ) on every compact set to another solution of this same equation. 

This follows from the boundedness of the sequence u o Tj and standard 
elliptic estimates. 

Lemma 3.2 Let L{z) be a bounded continuous function on HP which satis- 
fies L > c > and has the property that for any sequence of Mobius transfor- 
mations Tj, some subsequence of I°Tj converges uniformly on compact sets. 
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Let X be a vector field on HP with the analogous property, that for any such 
sequence Tj, the sequence of vector fields Xj obtained by pushing forward X 
by the action of Tj also has a convergent subsequence. Suppose that u is a 
solution to Au + X ■ Vu = I(z)u which is bounded below. Then inf u > 
and this infimum equals if and only if u = 0. The same conclusion holds 
if u is only a supersolution to this equation, but then we must also assume 
that uo Tj has a convergent subsequence for every sequence Tj. 

Proof: Let Zj be a sequence in HP such that u(zj) — ► infu. Choose a 
sequence of Mobius transformations Tj so that Tj (0) = Zj , and consider the 
sequence Uj = u o Tj, which solves the equation Auj + Xj • Vuj = IjUj, 
where Ij and Xj are the pullbacks of / and X by Tj. Using the hypothesis, 
choose a subsequence so that Ij converges locally uniformly to a continuous 
nonnegative function I. In some fixed neighbourhood of 0, Uj is also bounded 
above, and by elliptic estimates we may assume that Uj also converges in 
this neighbourhood to some function u. Now Au = Iu, I > c > and u 
attains a negative minimum at 0, which is impossible. □ 

3.1 The elliptic case 

The first symmetry result is proved by a 'stationary plane' reflection ar- 
gument, not too different from the standard moving plane method in the 
Euclidean case [5]. (We note, by the way, that [9] used moving planes to 
prove a symmetry result on certain bounded domains in hyperbolic space.) 

Theorem 3.3 Let u be a global solution of M.l\) on HP . 

i) If u is bounded and lim z _> 0O u{z) = rj, where either r] = +1 or n = —1, 
then u = n. 

ii) Suppose that f is Lipschitz, with Lipschitz constant L < (n — l) 2 /4. 
Suppose that \u(z)\/(l — \ z\) a < C for some a G (/?_, |(n — 1)). Then 
u = 0. The same conclusion is valid if f e C 1,s for some 5 > and 
a E (a-, \{n — 1)). Here a_ and (3- are defined in J l-4\ )- 

Proof: We prove ii) first. Consider u in the ball model. Fix any hyperplane 
P C HP and let R : HP -> HP be the reflection across P. Now define 
ur = uoR, which we consider as a function in one of the two components £1 
of HP \ P. Define wr = ur — u. We claim that wr > in 17. Interchanging 
the roles of the two components, we see that if this claim holds, then u is 
invariant with respect to reflections about P, and since P is arbitrary, we 
may conclude that u is constant. 
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Using that / is Lipschitz, we calculate that 



Awfl = /Cur) - f{u) < L\w R \. 

Suppose that wr < in some subdomain ST C Q. Then \wr\ = —wr there, 
hence Awr + Lwr < in Q'. Consider the function v = WR/4> a , where 4> a 
is the spherical function defined at the beginning of this section and a is as 
in the statement of this theorem. Then 

Vw R W R V<f> a 

VV = ; -7; , 



a • 



Av _ Aw R _ ^ Vwr • V4> a _ w R A(j) a 2 w R \V4> a \ 2 

so, using the equations for wt and 4> a , 

Av + X ■ Vv + (a(a - (n - 1)) + L)v < 0, 

where X = 2V(/> a /0 a . It is not hard to check, using the regularity of 
that X satisfies the hypothesis in Lemma [3.2i If a > where /?_ is defined 
in (II, 4p . then a (a — (n — 1)) + L < 0. Hence by this Lemma, v > in 17', 
so wr is also nonnegative, which is what we claimed. 

To prove the final assertion of this theorem, it suffices to assume that 
a > a_ if / is not just Lipschitz, but C 1,<5 for some 5 > 0, we invoke a 
regularity theorem from [TT], as implemented in [10]. To set this up, write 
(jl.ip as a perturbation of the solution which is identically 0: 

Au- f'(0)u = Q(u). 

Let A k ' u (H n ) denote the Holder space of order k + v, k € N, < v < 1, 
where all derivatives and difference quotients are computed with respect to 
the hyperbolic metric. Denote by p s A k,u the set of all functions u = p s v 
where v £ K k,u ; for convenience here and below, we denote by p = 1 — \z\ 2 
a smooth defining function for SqoHP, and we shall also assume that v = 5 
to reduce the number of different indices. The precise result we quote from 
these papers is as follows: fix A < (n — l) 2 /4, and let a± be the two indicial 
roots for A + A computed as in (jl.3p and (jl.4p ; then for any a S (a_,oq_), 
the mapping 

A + A : p a A 2 ' u (B. n ) — ► p a A°' u (M. n ) 

is an isomorphism. To apply this, observe first that since / G C 1,v , if \u\ < 
C(l - |z|) s for any s > 0, then |Q(«)| < C(l - |z|)( 1+l/ ) s . Now if |n| < Cp s 
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for s > a- and Ait = /(it), then it follows from scale-invariant Schauder 
estimates (namely, just ordinary Schauder estimates in balls of unit size 
with respect to hyperbolic distance) that it G p s A 3,,y (M n ) (i.e. precisely two 
derivatives more than the regularity of /). Now suppose u £ p a A 2,u is a 
solution; then Q(u) £ p a \ l + u ) h?' v , hence if a' := a(l + u) < a + , we deduce 
that u £ p a A 2,u . Continuing in this way a finite number of times implies 
eventually that u £ p a A 2,u for some a" > and now we may use the 
previous argument. 

To prove i) we must argue in a slightly different way since we do not 
wish to use any information about the size of the Lipschitz constant for /. 
To be definite, fix 77 = 1; the argument for 77 = — 1 is identical. 

First fix any geodesic 7 and let Pt denote the family of totally geodesic 
hyperplanes perpendicular to 7. Let R t be the reflection across Pt, and set 
itj = it o R t . Let denote the 'forward' component of EI n \ Pt, i.e. the 
component swept out by the P s , s > t, and the other component. We 
claim that when t <C 0, ut > u in Qf. To see this, note that on the one 
hand, for any given e>0, 1>u>1 — e when t is sufficiently negative. 
Suppose now that Ut < u in some subdomain $7' C Qf. Thus u(z) also lies 
in (1 - e, 1) for z £ fi'. 

The difference wt = ut — u satisfies the equation 



By assumption, f(s) is strictly increasing in some neighbourhood (1 — e, 1), 
and the argument of /' lies in this interval when z £ il', so I(z) > for 
z G 0,', and hence by Lemma 13.21 wt > 0. 

Now let to be the supremum of values in R for which itj > u in Qf. 
Clearly, unless u is constant, to must be finite since when i > 0, the roles 
of the components and ut and it are interchanged and we would have 
the reverse inequality. Either Pt must be a plane of symmetry for it or 
else ut Q > u in Qf Q . To see that this latter case is impossible, we must show 
that we can push Pt a bit further while still preserving this inequality, which 
contradicts that to is the supremum. For this it suffices to show that the 
normal derivative of wt is bounded below by a positive constant along all 
of Pt . If that were to fail, there would be a sequence of points qj £ Pt 
such that d v wt Q \ 0. Choose a sequence of Mobius transformations which 
preserve Pt and which move qj to a fixed point q. Applying Lemmas 13. II and 
13.21 we obtain a limiting function it; which satisfies Ait; = Iw with it; > 
and c^it; = at q, which violates the Hopf boundary point lemma. 




I(z) 
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We have now shown that one of the hyperplanes Pt orthogonal to the 
geodesic 7 must be a plane of symmetry for u. This is true for every geodesic 
7, and we claim that this implies that u is in fact identically constant. 

The proof of this last fact is actually slightly easier than in Euclidean 
space since in some sense there are 'more hyperplanes' in hyperbolic space. 
As shown above, for every (unoriented) geodesic 7 in HP, the solution u 
is symmetric around some hyperplane P 7 orthogonal to 7. Suppose that 
two of these hyperplanes, Pi and P2, do not intersect at any finite point 
of HP or on the asymptotic boundary S" 1-1 . Let c be the unique geodesic 
perpendicular to both of them. The composition of the reflections across 
Pi and P2 is a translation h along c which has the property that for any 
z £ HP, h k (z) — ► 00. However, u(h{z)) = u(z) for all z and u(h k (z)) — > 1, so 
u = 1. The same argument may be applied if P\ and P2 intersect at infinity, 
but not at an interior point; then h is a parabolic element, but it still has 
the property that every orbit {h k (z)} is unbounded, so the conclusion is the 
same. In fact, we see that the crucial property we need is that the group 
of motions T of HP generated by all the reflections P 7 across the various 
planes of symmetry P 7 satisfies that the orbit of any point z G M n under T 
is unbounded. Elementary hyperbolic geometry shows that this property is 
true for the group generated by any three reflections Pi, R2, R3, provided 
the corresponding planes do not all pass through the same point. In other 
words, either we can find some triplet of symmetry planes Pi, P2, P3 which 
do not pass through the same point, in which case we conclude that u = 1, or 
else all symmetry planes pass through some fixed point G HP. In this last 
case, we conclude that u is radial around this point. However, we showed in 
§2 that there are no nonconstant radial solutions to this equation, so once 
again, u must be constant. 

This concludes the proof in all cases. □ 

3.2 The parabolic case 

We have already noted that the only situation in which it is possible to 
expect a one-dimensional symmetry result when the boundary values are H p - 
invariant is when u — ► on 5' n ~ 1 \{(7} and u — > ±1 as z tends nontangentially 
to q. 

Theorem 3.4 Suppose that f is Lipschitz, with Lipschitz constant L < 
in — l) 2 /4. Let a > a-, and suppose that u is a global bounded solution of 
PUP satisfying \u(z)\/cr a < C on all o/HP, where the pole o/HP is placed 
at some point q £ 5 n ~ 1 , and moreover that u(z) — > 1 as z — > q in any 
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nontangential approach region in the ball. Then u is H p invariant, hence in 
the upper half-space model where q is at infinity, u(x,y) = U(x) where U is 
the function obtained in Proposition \2.2\ 

Proof: The proof is very close to that for part ii) of Theorem l3.31 Consider 
u in the upper half-space model, with the exceptional point q placed at 
infinity. Fix any vertical plane P given by i ■ y = c and let R denote the 
reflection across this plane and ur = u o R. As before, we set wr = ur — u, 
so that by hypothesis Awr = f{uR) — f(u) < L\wr\. Let ± be the two 
components of M n \ P. We claim that wr > in say. 

The proof is very similar indeed to what we had done before. If wr < 
in some subdomain $7', then we introduce v = WR/a a ; by the hypothesis, 
M < C, and v = on Q' n P, and v > — C in the entire subdomain. Since 
Av + X ■ Vi> < where X is a left-invariant vector field, we conclude that 
v > 0, a contradiction. Applying the same argument in the domain Q~ 
shows that ur = u, so P is a plane of symmetry. Since every vertical plane 
is now a plane of symmetry, u is independent of y, and hence a function of 
x alone. 

It is very important here that wr > —Ca a , so that we may control the 
behaviour of v as (x, y) — > oo nontangentially, e.g. as y — > oo with x fixed. 
□. 

Note that even for ODE solutions, tangential limits of u at q can take 
on other values in (—1,1). 

3.3 Hyperbolic case 

The corresponding symmetry result in the hyperbolic case is less standard; 
the proof uses what could be called a 'stretching method', in analogy to the 
'sliding method' in the Euclidean setting 

Theorem 3.5 Let u be a solution of U.l\) satisfying the asymptotic bound- 
ary condition that u(z) converges to ±1 at any point in the interior of the 
two hemispheres S C cffl™. Let HP -1 be the totally geodesic subspace with 
boundary S- n S+, and t the signed distance function to this subspace in M n . 
Suppose that f is a nonnegative C 1 function onM. such that /(l) = /(— 1) = 
and /'(l) > 0, /'(— 1) > 0. Thenu(z) depends only ont, i.e. u(z) = U(t(z)), 
where U satisfies the ordinary differential equation 

U"(t) + (n - l)^-U'(t) = f(U), hni U(t) = ±1. 
cosh t t^±oo 
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Proof: Let P_ G S- and P+ G S+ be arbitrary points, and suppose that 7(t) 
is any curve of constant geodesic curvature in HP such that lim£_ > -i- 00 7(i) = 
P±. We shall first prove that u o j(t) is (weakly) monotone increasing in t. 
The conclusion is not difficult to prove from this. 

For this step, we work in the upper half-space model. Compose with 
a Mobius transformation so that P_ = and P+ = oo; the images of S- 
and S+ (which we denote by the same symbols) are then some ball in IP -1 
containing 0, but not necessarily at its center, and the exterior of this ball 
(union oo), respectively, and the curve 7 is transformed into some ray in the 
upper half-space emanating from 0. 

Fix 5 > such that f(s) > when s G [-1, -1 + S) U (1 - 6, 1], and 
choose A G (0, 1) so that u(z) < -1 + 5 for z G D-{A) := {z G : \z\ < A} 
and u{z) > 1 - 5 for z G D+(A) := {z G R™ : \z\ > l/A}. We let diD±(A) 
denote that portion of dD±(A) lying in HP. 

Now choose any R > l/A 2 and define ur(z) = u{Rz). We claim that 
ur{z) > u(z) for all z G M n . We verify this in D-(A) and M n \ D-(A) 
separately. First note that if z G dD-(A), then nij(z)>l — <5>— 1 + 5 > 
u(z). Next, suppose that the subset Q = {z G D_(vl) : Ur(z) < u(z)} C 
D_(j4) is nonempty, and define wr = ur — u. Then in f2, 

Awr = I(z)wr, I(z)= [ f'(u(z)+t( U R(z)-u(z)))dt, (3.11) 

Jo 

while wr = on the entire boundary of O, including the portion lying in 
dW n . Furthermore, if z G f2, then ur(z) < u(z) < —1 + 5, and hence 
u{z) + t(u fl (*) - < -1 + 5 for all < t < 1. Since 5 < 1, !(*) > 

on Q. We can now apply the maximum principle to (|3.1ip to conclude that 
ur — it > in f2, which is a contradiction to the definition of this set, and 
hence 17 = 0. Finally, if z G M n \ D_(A), then ur(z) > 1 — 5, and a similar 
argument proves that ur> u holds outside of D_{A) as well. This proves 
the claim. 

Now define R = ml{R > 1 : u R (z) > u(z) Vz G HP}. We have just 
established that this set is nonempty, so that R G [1, 00) is well defined. We 
claim that R = 1. 

Suppose by way of contradiction that R > 1. Define u = Ur, wr = ur—u 
and w = u — u. By continuity, w > on HP, and w satisfies an elliptic 
equation, as in (|3.1ip . with coefficient I(z). Decomposing I = J + — we 
see that Aw = (/+ — lJ)w < I+w, so by the strict minimum principle for 
supersolutions, w > on H n . 

If R > 1, there would exist a sequence Rk R and points 2^ G H n so 
that the function WR k := wt attains a strictly negative absolute minimum 
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value at Zk- If some subsequence of the points Zk were to converge to a point 
z 6 HP, then we would have w{z) = 0, which we have shown cannot happen. 
Thus Zk diverges in HP. 

To renormalize, choose a sequence of Mobius transformations, T^, so that 
Tk(zk) is equal to some fixed point z £ HI™. Denoting wj. o Tk by w^, then 
as before, each Wk satisfies an equation Awk = Ik(z)wk, where Ik = Ik° Tk', 
furthermore, inf Wk = Wk(z) = —tk S 0. 

Using Lemmas 13.11 and 13.21 again, we may assume that % — > u, w/- w 
and Ik —* I, and 

Aw = Iw, w > 0, w(z) = 0. 

Arguing as above, we may apply the strict minimum principle for superso- 
lutions again to conclude that w = 0. This shows that I = f'(u); however, 
since this function does not necesssarily have a sign, we must proceed further. 
In fact, we derive a contradiction by considering the function Wk = ^Wk 
(where — e& = infu^). Clearly AWk = IkWk, but now min Wk = — 1. As 
before, we can extract a subsequence converging in C 2,a on every compact 
set to some function W, where 

AW = IW, MW = W{z) = -1. (3.12) 

Now let us show why such a limit cannot exist. There are two cases to 
consider, depending on the limiting behaviour of the sequence Zk in HP\HP: 

i) Suppose that Zk converges to a point Zoq which lies in either of the 
open half-balls S± where u is continuous and equals ±1. In this case, 
I(z) = /'(limtifc(z)) = /'(±1) > 0, so (I3T21) is impossible. 

ii) Suppose on the other hand that Zk converges to some point in the 
intersection of the two hemispheres S+ n S-. Now recall that (using 
upper half-space notation again), since Rk > c > 1 for all k, then 
RkZk must converge to a point in S+, so u(RkZk) — > 1. However, by 
definition, u(RkZk) — u(zk) < 0, so u(zk) — ► 1 as well. (This is not 
incompatible with z^ lying on 5+ fl S- , but rather, indicates that it 
approaches this intersection tangent to the S+ side.) Thus here too 
we conclude that I(z) > 0, at least in a neighbourhood of z, so f|3. 12|) 
is again impossible. 

We have now ruled out the possibility that R > 1, so R = 1. In other 
words, we have proved that ur > u for all R > 1; hence 

u(Rz) = rd r u(z) > 
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for every z G HP; here r is the Euclidean radial upper-half-space variable. 
This shows that u is nondecreasing along any ray emanating from the origin, 
i.e. along any curve of constant geodesic curvature joining P_ to P+. 

We now prove that u is a function of only one variable. Working again in 
the ball model, suppose that S_ and S + are the lower and upper hemispheres 
of the boundary, respectively. Let P be any point in the interior of the ball, 
and let E be the spherical cap which passes through P and S- H S+. Thus E 
is a hypersurface equidistant from the totally geodesic copy of HP -1 which 
has the same asymptotic boundary S + n S-. Let II be any two-dimensional 
plane passing through the origin of the ball and the point P. Then LT n E 
is a curve of constant geodesic curvature 7 in HP passing through P and 
limiting on two points Q, Q' G S- Pi S + . It is easy to see geometrically (for 
example by tilting E slightly, but still ensuring that it passes through P), 
that we can approximate 7 by two sequences of curves of constant geodesic 
curvature, 7j~(i) and 7^~(t) such that 7^(0) = P for all j, 

lim 7- (t ) lim 7 + (t) G SL , lim 7 7 (f ) lim 7+ (i)eS + , 

r^— 00 J t— >+oo J t— >+cxd J t— >— oo J 

and finally 

TpE 9 X = lim ( 7 7)'(0) = - lim (7, + )'(0). 

Since uo^j(t) and uo^(t) are both nondecreasing, we see that Vup-X = 0. 
However, X can be chosen arbitrarily in TpE, which shows that Vu(P) is 
orthogonal to S. 

We have now proved that if E^ is the foliation of HP by hypersurfaces 
which are of (signed) distance t from the totally geodesic copy of HP -1 with 
boundary S- H S + , then each E t is a level set of u. In other words, u is a 
function of the distance t alone. This concludes the entire proof. □ 

4 Non- symmetric solutions 

We have already indicated that the decay hypotheses in parts i) and ii) of 
Theorem 1.2 are necessary since in fact there is a very large class of non- 
symmetric solutions of (II. ID the elements of which satisfy the bounds 

a) \u(z)\ < C(l - \z\ 2 ) a ~, 

b) \u(z)\ <Cx a ~, and lirn^oo u(x, y) = 1 

respectively. In this section we give a brief sketch of the construction of 
these solutions under the slightly stronger assumption that / G C 1,<5 (this is 
probably not necessary but is used in the proof). 
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This construction is a simple perturbation argument using the implicit 
function theorem. Suppose that no is a known solution of (jl.ip . Then any 
nearby solution u may be decomposed as u = no + v where v is small, both 
globally and perhaps also in terms of its asymptotic decay rate. Expanding 
/ into its first order Taylor polynomial and an error term, the equation 
An = f(u) can be rewritten as 

Av = f(uo)v + Q(u ;v), 

where Q at least has the property that if \v\ < e then |Q(w)| < Ce 1+S , and 
the notation v i— ► Q{uq;v) is meant to indicate the dependence of Q on uq. 

We shall actually use a slightly finer decomposition. Let <f> any element 
of the nullspace of A — f'(uo); in the examples below, <p will have the same 
asymptotic behaviour as no, but will not decay faster. Now write u = 
uo + <j) + w, where w lies in some space of functions which decay faster than 
uo + (j) at infinity. We rewrite the equation one last time as 

(A - f{ua))w = Q(u ; (/> + w). (4.13) 

The existence of nearby solutions can therefore be deduced from the 
mapping properties of the operator L = A — f'(uo) and from the asymptotic 
properties of solutions of Lcf) = 0. We now describe these in the two cases of 
interest. The reason that this sort of analysis works is because we consider 
perturbations of solutions no which are asymptotically equal to on most 
or all of OoqW 1 and because is an unstable critical point of the potential 
function F. We shall use the weighted Holder spaces p s A k ' u introduced in 
§3.1. 

For the elliptic case, let no = 0. Recall the Poisson transform P on 
hyperbolic space. It is well-known that P provides a correspondence between 
C 2 ^{S n - 1 ) and the space of solutions of (A - f (0))<£ = with <j) ~ <j)op a ' , 
foeC^iS 11 - 1 ). 

Proposition 4.1 Let no = and fix a £ (a_,a+) ; where a± are described 
in M.3\) and {l-4\) - There is a neighbourhood of U C C 2 ' /i (S' n_1 ) and a 
map G : U -> V C p a k^ v {W) so that all solutions u G p a - k 2 > u to [JJ^ j 
close to are of the form u = P(<f>o) + G(4>o)- 

Proof: First note that (4>o,w) t— > Q(O;P(0o) + n) is a continuous mapping 
from U x V to p a A 2 ' u (H. n ) provided a < (1 + 5)a_. (Here U and V are small 
balls in their respective function spaces, as in the statement of this result.) 
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Furthermore, as described in §3.1 and proved in [TT], since A = —f'[0) < 
(n-l) 2 /4, 

A + A : p a A 2 > u (U n ) — > p a A°' u (M n ) 
is an isomorphism. We now rewrite (|4.13|) as 

w = (A + A)- 1 Q(O;P(0 O )+^) 

and, for each 0o sufficiently close to 0, obtain a solution by showing that 
the map on the right is a contraction mapping. □. 

In the parabolic case we follow a similar strategy. Since we are not 
trying to obtain the most general possible result, but only wish to exhibit 
the abundance of possible solutions, we simplify the analysis by restricting 
the allowable perturbations to have a very special form. Let uo(x) denote 
the ODE solution obtained in Theorem 1.1, part ii), and consider nearby 
solutions u(x, y) which are invariant with respect to the integer lattice T = 
Z n_1 acting on y E M n_1 . In other words, we consider the problem on the 
quotient M = HP/T. This manifold M is diffeomorphic to the product 
K + xT,T = T n ~ l = W 1 ' 1 /!^- 1 - it has two ends, one a 'large' conformally 
compact end to which the same sort of analysis as above applies, and the 
other a finite volume hyperbolic cusp end. As a further simplification, we 
shall consider perturbations of the form u(x,y) = uq(x) + v(x,y) where v 
is not only periodic in y, hence descends to this quotient, but in addition, 
fjv(x,y) dy = for all x > 0. (This is sometimes called the cusp form 
condition.) 

Let 0o £ C 2,a (T) be any function with vanishing integral. There is 
once again a Poisson transform P which maps 4>q to the solution <j>(x, y) 
to (A — f'(uo))4> = with (j>(x,y) ~ <po(y)x a - as x \ 0. This solution </> 
satisfies the cusp form condition, and decays exponentially as x f oo. All of 
these facts may be proved by separating variables and using standard ODE 
arguments. 

The next step is to construct a parametrix G for L := A — f'(uo) by 
pasting together parametrices for A — /'(0) near x = and for A — /'(l) for 
x ^> 0. Using the same sort of analysis from [TT] as employed in the previous 
result, there is an inverse Go for A — /'(0) on all of M; similarly, there is 
also an inverse Goo for A — f'(l) since f'(l) > 0. Now choose a partition 
of unity xo + Xoo = 1 where xo = for x 3> and Xoo = for x < 1, say. 
Also choose functions xo> and Xoo such that Xj = 1 on the support of Xji 
j = 0, oo, and xo vanishes for x ~S> 0, Xoo = for x < 1/2, say. Then 

G = XoGqXo + 
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is a parametrix for L in the sense that LG = I — R where R is a smoothing 
operator which induces extra decay. We can choose all of these operators 
to preserve the subspace of functions on M which satisfy the cusp-form 
condition. 

To be more precise about mapping properties, set (x) = y/1 + x 2 , and 
consider weighted spaces of the form (x/(x)) a (x)^C k ' u (M). These functions 
grow or decay like x a and x 13 as x \ 0, x /* oo, respectively. Then 

G : (x/(x)) a (xfA°' a — ► (x/(x)) a (x)^A 2 ' a 

is bounded for a £ (a_, a+) and (3 £ (—y/f'(l), y /'(l)). Restricting to the 
subset of functions satisfying the cusp-form condition, we could even let (3 
be an arbitrary real number. Furthermore, 

R : (x/(x)r(xfA° > a — (x/ixyr^ixf-'A^, 

from which its compactness follows easily. All of this shows that 

L : (x/(x)r(x)PA 2 > a — » (x/(x))«(x)^ Q , 

is Frcdholm. It remains to show that this map is actually an isomorphism, 
which is true and can again be verified using separation of variables and 
ODE arguments. 

The rest of the argument is exactly the same as in the previous case. 
This proves the 

Proposition 4.2 There is an infinite dimensional space of solutions of 
Au = f(u) on HP /r which are small perturbations of the ODE solution 
uo(x). Each such solution has an asymptotic expansion of the form u(x, y) ~ 
(c+(po(y))x a - as x \ 0, for all sufficiently small <fio satisfying the cusp-form 
condition. 
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